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Abstract. We study the number of limit cycles and the bifurcation diagram in the 
Poincare sphere of a one-parameter family of planar differential equations of degree 
five X = Xb(x) which has been already considered in previous papers. We prove 
that there is a value b* > such that the limit cycle exists only when 6 G (0,6*) and 
that it is unique and hyperbolic by using a rational Dulac function. Moreover we 
provide an interval of length 27/1000 where b* lies. As far as we know the tools used 
to determine this interval are new and are based on the construction of algebraic 
curves without contact for the flow of the differential equation. These curves are 
obtained using analytic information about the separatrices of the infinite critical 
points of the vector field. To prove that the Bendixson-Dulac Theorem works 
we develop a method for studying whether one-parameter families of polynomials 
in two variables do not vanish based on the computation of the so called double 
discriminant. 
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1. Introduction and main results 
Consider the one-parameter family of quintic differential systems 
( X = y, 

\y = -X + {a - x'^){y + y^), a G M. ^' 

Notice that without the term y^, coincides with the famous van der Pol system. 
This family was studied in [21] and the authors concluded that it has only two 
bifurcation values, and a*, and exactly four different global phase portraits on the 
Poincare disc. Moreover, they concluded that there exists a* G (0, i^Ott^/IG) ~ 
(0, 1.77), such that the system has limit cycles only when < a < a* and then if 
the limit cycle exists, is unique and hyperbolic. Later, it was pointed out in [11] 
that the proof of the uniqueness of the limit cycle had a gap and a new proof was 
presented. 

System ([1]) has no periodic orbits when a < because in this case the function 
x^ + y"^ is a global Lyapunov function. Thus, from now on, we restrict our attention 
to the case a > and for convenience we write a = b"^, with b > 0. That is, we 
consider the system 

{ y = -x+ - x^){y + y^), 6 G M+ U {0}. 

Therefore the above family has limit cycles if and only if 6 G (0, b*) with b* = 
^/a* and b* G (0, ^97rVl6) ^ (0, 1.33). Following [21] we also know that the value 
b = corresponds to a Hopf bifurcation and the value b* to the disappearance of 
the limit cycle in an unbounded polycycle. By using numerical methods it is not 
difficult to approach the value b*. Nevertheless, as far as we know there are no 
analytical tools to obtain the value b*. This is the main goal of this paper. 

We have succeed in finding an interval of length 0.027 containing b* and during 
our study we have also realized that there was a bifurcation value missed in the 
previous studies. Our main result is: 

Theorem 1.1. Consider system ([2]). Then there exist two positive numbers b and 
b* such that: 

(a) It has a limit cycle if and only if < b < b* . Moreover, when it exists, it 
is unique, hyperbolic and stable. 

(b) The only bifurcation values of the system are 0,b and b* . In consequence 
there are exactly six different global phase portraits on the Poincare disc, 
which are the ones showed in FigureUi 

(c) It holds that 0.79 < b < b* < 0.817. 

The phase portraits missed in [21] are (ii) and (iii) of Figure [1] 
The key steps in our proof of Theorem 11.11 are the following: 

• Give analytic asymptotic expansions of the separatrices of the critical 
points at infinity, see Section [21 

• Use these expansions to construct explicit piecewise rational curves, and 
prove that they are without contact for the fiow given by These 
curves allow to control the global relative positions of the separatrices of 
the infinite critical points, see Section O 
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(o)a<0 (i)0<b<b (ii)h = h 




(iii)b<b<b* (iv) b = b* (v) b > b* 



Figure 1. Phase portraits of systems ([T]) and ([2]). When a > 0, 
then b = ^Ja. 

• Provide an alternative proof of the uniqueness and hyperbohcity of the 
hmit cycle, which is based in the construction of an explicit rational Dulac 
function, see Section |H 

By solving numerically the differential equations we can approach the bifur- 
cation values given in the theorem, see Remark 12.61 We have obtained that 
b ^ 0.8058459066, b* ^ 0.8062901027 and then b* -b ^ 0.000444. As we have said 
the main goal of this paper is to get an analytic approach to the more relevant 
value b*, because it corresponds to the disappearance of the limit cycle. 

Although all our efforts have been focused on system ([2]), the tools that we 
introduce in this work can be applied to other families of polynomial vector fields 
and they can provide an analytic control of the bifurcation values for these families. 

As we will see, our approach is not totally algorithmic and following it we do 
not know how to improve the interval presented in Theorem 11.11 for the values b 
and b* . 

One of the main computational difficulties that we have found has been to prove 
that certain polynomials in x, y and b, with high degree, do not vanish on some 
given regions. To treat this question, in Appendix II we propose a general method 
that uses the so called double discriminant and that we believe that can be useful 
in other settings, see for instance [H [22] . In our context this discriminant turns 
out to be a huge polynomial in 6^ with rational coefficients. In particular we need 
to control, on a given interval with rational extremes, how many reals roots has 
a polynomial of degree 965, with enormous rational coefficients. Although Sturm 
algorithm theoretically works, in practical our computers can not deal with this 
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problem using it. Fortunately we can utilize a kind of bisection procedure based 
on the Descartes rule ([12]) to overcome this difficulty, see Appendix I. 

2. Structure at infinity 

As usual, for studying the behavior of the solutions at infinity of system ([2]) 
we use the Poincare compactification. That is, we will use the transformations 
{x,y) = {l/z,u/z) and {x,y) = {v/z,l/z), with a suitable change of time to 
transform system ([2]) into two new polynomial systems, one in the {u, z)-plane 
and another one in the (f , z)-plane respectively (see [2] for details). Then, for 
understanding the behavior of the solutions of near infinity we will study the 
structure of the critical points of the transformed systems which are localized on 
the line z = 0. Recall that these points are the critical points at infinity of system 
02]) and their separatrices play a key role for knowing the bifurcation diagram of 
the system. In fact, it follows from the works of Markus [16] and Newmann [T7] 
that it suffices to know the behavior of these separatrices, the type of finite critical 
points and the number and type of periodic orbits to know the phase portraits of 
the system. We obtain the following result: 




Figure 2. Separatrices at infinity for system (I2j). 

Theorem 2.1. System has six separatrices at infinity, which we denote by 
Si,S2, S3,S[, S'2 and S'^, see Figure\^ Moreover: 

{i) Each is the image of Sk under the transformation {x,y) (— x, —y). 
{ii) The separatrices S2 and S3 near infinity are contained in the curve {y — 
4>{x) = 0} where (f){x) = (f){x — b)/{x — 6)^, (f){u) is an analytic function at 
the origin that satisfies 

In particular, S2 corresponds to x '^b and S3 to x > b. 
[Hi) The separatrix Si near infinity is contained in the curve {y — (f{x) = 0} 
where (p{x) = ip{l/x) and (p is an analytic function at the origin that 
satisfies 

(p{u) = -u- (b^ - l)u^ - (b^ - 3b^ + 2)u^ + O {u') . (4) 

Remark 2.2. In the statements {ii) and {Hi) of Theorem \2.1\ the Taylor expan- 
sions of the functions cf) and (p can be obtained up to any given order. In fact, in 
Section\^we will use the approximation of (j) until order 16. 
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(iv) (v) 
Figure 3. Relative position of the separatrices of system ([2]). 

Corollary 2.3. All the possible relative positions of the separatrices of system ([2]) 
in the Poincare disc are given in Figure\^ 

To prove the above theorem we need some preliminary lemmas. 

Lemma 2.4. By using the transformation {x,y) = {l/z,u/z) and the change of 
time dt/dr = \jz^ system ([2]) is transformed into the system 

where the prime denotes the derivative respect to r. The origin is the unique 
critical point of ^ and it is a saddle. Moreover the stable manifold is the u-axis, 
the unstable manifold, Si, is locally contained in the curve {u — ip{z) = 0}, where 
ip{z) is an analytic function at the origin that satisfies 

^(z) = -z^ - {}? - l)z^ - {b* - 36' + 2)/ + 0(z«), (6) 

see Figure\^ 




Figure 4. Phase portrait of system 
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Proof. From the expression of ([5]) it is clear that the origin is its unique critical 
point. For determining its structure we will use the directional blow-up since the 
linear part of the system at this point vanishes identically. 

The u-directional blow-up is given by the transformation u = u, q = z/u; and 
by using the change of time dt/dr = v?, system ([5]) becomes 

u = —u — {1 — h'^u^)uq^ ~ (1 — h'^u)u'^q^ — u^q^, , . 

q = q+(l-b'^u^)q^+{l-b^u)uq^. 

This system has a unique critical point at origin and it is a saddle with eigenval- 
ues ibl. 

The z-directional blow-up is given by the transformation r = u/z, z = z. Doing 
the change of time dt/dr = —z"^, system ([5]) becomes 

r = z + {l-b'^z'^){r + r^), 

4 \°) 
z = rz . ^ ' 

This system has a unique critical point at the origin which is semi-hyperbolic. We 
will use the results of [21 Theorem 65] to determine its type. By applying the 
linear change of variables r = + i], z = C, system ([8]) is transformed into 

7] = 7]- N{^,r]), 

where N{^,7]) = (r/ - - + iv - O^i^'^^'^ - !)• It is easy to see that if 

V = "^{0 is the solution of 77— A^(^, 77) = passing for the origin, then n(^) = — (6^ — 
1)^^ - {b^ - + 2)^5 + 0(^^). Thus - = + 0{^'^). Therefore from 
[21 Theorem 65] we know that the origin is a semi-hyperbolic saddle. Moreover, 
its stable manifold is the ?7-axis and its unstable manifold is given by 

= -(^2 _ 1)^3 _ (^4 _ 3^2 ^ 2)^5 ^ Q^^Vy 

In the plane (r, z) the local expression of this manifold is 

r = -z-{b^- 1)^3 _ _ 3^2 ^ 2)^5 + O(z^). 

Finally, in the {u, z)-plane the unstable manifold is contained in the curve 
and from the analysis of phase portraits of systems ([7D and ([HD we obtain that the 
local phase portrait of system (jS]) is the one given in Figure HI □ 

Lemma 2.5. By using the transformation {x,y) = {v/z,l/z) and the change of 
time dt/dr = 1/z^ system is transformed into the system 



v' = v{l + z^)iv^-bh^) + {l+v^)z^ 
z' = z{l + z^)lv'^ - b'^z'^) + vz' 



where the prime denotes the derivative respect to r. System has a unique 
critical point at the origin and its local phase portrait is the one showed in Figure\^ 
Moreover, the separatrices S2 and S-s are locally contained in the curve {v — g{U) = 
0} where U = z/v — l/b and g{U) is an analytic function at the origin that satisfies 
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Figure 5. Topological local phase portrait of system All the 
solutions are tangent to the f-axis but for aesthetical reasons this 
fact is not showed in the figure. 

Proof. From the expression of system it is clear that the origin is its unique 
critical point. As in Lemma (2.41 we will use the directional blow-up technique to 
determine its structure since the linear part of the system at this point is identically 
zero. 

It is well-known, see [2], that since at the origin z'v — v'z = —z^ + 0{z^), all 
the solution, arriving or leaving the origin have to be tangent to z = 0. So it 
suffices to consider the w-directional blow-up given by the transformation v = v, 
s = z/v. Performing it, together with the change of time dt/dr = —v^, system 
(jnD is transformed into 

\ s = s^. ^ ^ 

This system has not critical points. However, by studying the vector field on the 
s-axis we will obtain relevant information for knowing the phase portrait of system 
(jni). If s = then v = —1 and s = 0, that is, the v axis is invariant. If = then 
V = — l + and s = s^, this implies that = if s = ±1/6. In addition, a simple 
computation shows that i) > at the points (0, ±1/6). Therefore the solutions 
through these points are as it is showed in Figure ISi(a), and by the continuity 
of solutions with respect to initial conditions, we have that the phase portrait of 
system ([9]), close to these points, is as it is showed in Figure El (b). 




(a) (b) 

Figure 6. Local phase portrait of system f lTT]) . 

Then by using the transformation [v, z) = {v, sv) and the phase portrait showed 
in Figure [ni(b) we can obtain the phase portrait of system Recall that the 
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mapping swaps the second and the third quadrants in the w-directional blow-up. 
In addition, taking into account the change of time dt/dr = —v^ it follows that 
the vector field in the first and fourth quadrant of the plane {v, z) has the opposite 
direction to the showed in the (w, s)-plane. Therefore the local phase portrait of 
© is the showed in Figure [51 




Figure 7. Transformation between system (ITT]) and system ([H]). 

To show that the separatrices S2 and 1S3 are contained in the curve (|T0|) we 
proceed as follows. First, we will obtain the curve that contains the solution 
through the point (0, 1/6) in the plane (w, s). Second, by using the transformation 
(f , z) = (f , sv) we will obtain the corresponding curve in the {v, 2;)-plane and we 
will show that such curve is exactly the curve given by ffTOl) . 

Since s is positive in (0, 00), the solution through the point (0, 1/6) (respectively 
(0,-1/6)) is contained in the curve {v — g{s) = 0} (respectively {v — g{s) = 
0}), where g{s) (respectively g{s)) is an analytical function defined in an open 
neighborhood of the point, moreover it is clear that g{l/b) = and g'{l/b) = 0. 
Consider the Taylor series of g{s) around (1/6): 

i=2 

Since the curve {v — g{s) = 0} is invariant then (V(w — g{s)),X) = at all the 
points of {v — g{s) = 0}, where X is the vector field associated to system ffTTl) . 
Thus, we have a function, {V{v — g{s)), X), for which all its coefficients have to be 
zero. From this observation we obtain linear recurrent equations in the coefficients, 
g^^\l/b) of g{s). Simple computations show that the first 3 terms of the Taylor 
series of g{s) are: 

Thus, in the plane {v, z), the curve corresponding to {v — g{s) = 0} is 

w (f - i)' - f (f ((f -ir)=o}. 

Finally, if f/ = ^/t; - 1/6, we obtain (do]). □ 

Remark 2.6. The proof of the above lemma gives a natural way for finding a 
numerical approximation of the value 6*. Notice that in the coordinates {v, s) the 
point (0, 1/6) corresponds to both separatrices S2 and S3. Since it is a regular point 
we can start our numerical method (we use a Taylor method) without initial errors 
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and then follow the flow of the system, both forward and backward for given fixed 
times, say t~^ > and t~ < 0. We arrive to the points [v"^, s^) with ^ for 
t = t^, respectively. These two points have associated two different points (x^, y^) 
in the plane {x,y), because of the transformation {v,s) = {x/y,l/x). Now, we 
integrate numerically the system ([2]) with initial conditions {x^,y^) to continue 
obtaining approximations of the separatrices S2 and S3, respectively. The next step 
is to compare the points of intersection x~^ = x~^{b) < and x~ = x~{b) > of 
these approximations with the x-axis. 

We consider the function b — )■ 11(6) := x~^{b) + x~{b) and we use the bisection 
method to find one approximate zero of U. Note that if 11(b) = then S2 = 
S3 and by the symmetry of the system S^ = S2, and therefore b* = b. Taking 
bo = 0.8062901027, t+ = 0.05 and r = -0.5 we obtain that x+{bo) + x^{bo) ^ 
-4.58036036 x lO-^^ and so b* ^bo. 

Following a similar procedure, but now using Lemma \2.4\ to have an initial 
condition almost on Si, we get that b ~ 0.8058459066. 

Proof of Theorem \2.1[ (i) The result follows because system (j2]) is invariant by the 
transformation {x,y) — )■ {—x,—y). 

(ii). From (|T0|) and by using the change of variables {v,z) = {x/y,l/y) we 
obtain that the separatrices ^2 and S3 are contained in the curve 



\y \x bj 3 \x bj 3 \x b 
or equivalently 




{y - (f){x) = 0} , (13) 



where 

0(x 



bHi-iY-fbni-ir+fb^{i-ir+o[{i-iy 

We can write the function (t>(x) as 



where 

b'^x^ 

The function 0i(x) is analytical ai x = b and it is not difficult to see that it has 
the following Taylor expansion 

- , , 1 (x-b) [x-bf 359(x - bf 



—^-6 35^ ' 963 276^ ,0((.-6)). 

Then 0141) can be written as 

?w = t,;^ - + 9^ - - ') + - i-n. 

Hence from (IT^ and taking 0(x) = 0(x — b)/{x — b)"^ we complete the proof. 

The proof of (iii) follows by applying the previous ideas, considering the expres- 
sion given by (|6]) and the change of variables [u, z) = [y /x,l/x). □ 
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3. Proof of Theorem 11.11 

We start proving a preliminary result that is a consequence of some general 
properties of semi-complete family of rotated vector fields with respect one pa- 
rameter, SCFRVF for short, see [71 [TB]. 

Proposition 3.1. Consider system ([2j) and assume that for b = b > it has no 
limit cycles. Then there exists < b* < b such that the system has limit cycles 
if and only if b & {0,b*). Moreover, for b = b* its phase portrait is like (iv) in 
Theorem \l.l\ and when b > b* it is like (v) in Theorem \l.l\ 

Proof. It is easy to see that the system has a limit cycle for 6 > 0, which appears 
from the origin through an Andronov-Hopf bifurcation. 

If we denote by Xb{x,y) = {Pb{x,y),Qb{x,y)) the vector field associated to ([2]) 
then 

A,,, tan (9i!(^\ = P^i^,yf-^-Q^i^,yf-^ 

962^'^ \Pb{x,y)J P^{x,y) + Ql{x,y) 

Pi{x,y) + Ql{x,y) " 

This means that system ([2]) is a SCFRVF with respect to the parameter b^. 

We will recall two properties of SCFRVF. The first one is the so called non- 
intersection property. It asserts that if 71 and 72 are limit cycles corresponding to 
different values of b, then 71 fl 72 = 0. 

The second one is called planar termination principle: [121 [20] if varying the 
parameter we follow with continuity a limit cycle generated from a critical point 
p, we get that the union of all the limit cycles covers a 1-connected open set U, 
whose boundaries are p and a cycle of separatrices of Xi,. The corners of this cycle 
of separatrices are finite or infinite critical points of Xh- Since in our case Xb only 
has the origin as a finite critical point we get that U has to be unbounded. Notice 
that in this definition, when a limit cycle goes to a semistable limit cycle then we 
continue the other limit cycle that has collided with it. This limit cycle has to 
exist, again by the properties of SCFRVF. 

If for some value of b = b > the system has no limit cycle it means that the 
limit cycle starting at the origin for b = 0, has disappeared for some b*, < b* <b 
covering the whole set U. Since U fills from the origin until infinity, from the non 
intersection property, the limit cycle cannot either exist for 6 > 6*, as we wanted 
to prove. 

Since for 6 > the origin is a repellor, by Corollary l2.3l we know by the Poincare- 
Bendixson Theorem that the phase portraits (i),(ii) and (iii) in Figure [1] have at 
least one limit cycle. Then, the phase portraits for b > b* have to be like (iv) or 
(v) in the same figure. Since the phase portrait (iv) is the only one having a cycle 
of separatrices it corresponds to b = b*. Again by the properties of SCFRVF, the 
phase portrait (iv) does not appear again for b > b*. Hence, for b > b* the phase 
portrait has to be like (v) and the proposition follows. □ 

Remark 3.2. In Lemma \4^ we will give a simple proof that when 6 = 1 system ([2]) 
has no limit cycles, based on the fact that for this value of the parameter it has the 
hyperbola + 1 = invariant by the flow. From the above proposition it follows 
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that b* < 1. This result already improves the upper bound of b* , given in [21] . 
^97r2/16 ^ 1.33. Theorem\r^ improves again this upper bound, but as we will 
see, the proof is much more involved. 

Proof of Theorem \l.l[ Recall that for a < the function V{x, y) = + y^ is a 
global Lyapunov function for system ([1]) and therefore the origin is global asymp- 
totically stable. Then it is easy to see that its phase portrait is like (o) in Figure [H 
To prove the theorem we list some of the key points that we will use and that 
will be proved in the forthcoming sections: 

(Ri) System ([2]) has at most one limit cycle for b G (0, 0.817] and when it exists 

it is hyperbolic and attractor, see Section HI 
(R2) System ([2]) has an odd number of limit cycles, with multiplicities taken 
into account, when b < 0.79 and the configuration of its separatrices is like 
(i) in Figure [31, see Proposition 15.11 
(R3) System ([2]) has an even number of limit cycles, with multiplicities taken 
into account, when b = 0.817 and the configuration of its separatrices is 
like (v) in Figure [3l see again Proposition 15.11 
The theorem for b > b* is a consequence of Proposition 13. 1[ Notice that again 
by this proposition and (R3), b* < 0.817. Hence, the limit cycles can exist only 
when b G (0, b*) C (0, 0.817] and by (Ri) when they exist then there is only one 
and it is hyperbolic and attractor. 

As a consequence of (R2) and the uniqueness and hyperbolicity of the limit 
cycle we have that the phase portrait for b < 0.79 is like (i) in Figure [H 

To study the phase portraits for the remaining values of b, that is 6 G (0.79, b*), 
first notice that all of them have exactly one limit cycle, which is hyperbolic 
and stable. So it only remains to know the behavior of the infinite separatrices. 
We denote by X2{b) and x'^^i^b) the points of intersection of the separatrices S2 
and 1S3 of system ([2]) with the x-axis (when they exist), see also the forthcoming 
Figure [13 Notice that for b > b*, x'^{b) < ^2(6) < and x'^{b*) = ^2(6*) < 0. 
The properties of the SCFRVF imply that X2{b) is monotonous increasing and 
that x'^{b) is monotonous decreasing. Hence for b ^ b* the phase portrait of the 
system is like (iii) in Figure [1] Since we already know that for 6 = 0.79 the phase 
portrait is like (i), it should exists at least one value, say b = b, with phase portrait 
(ii). Since for SCFRVF the solution for a given value of b, say b = b, becomes a 
curve without contact for the system when 6 7^ 6, we have that the phase portraits 
corresponding to heteroclinic orbits, that is (ii) and (iv) of Figure [H only appear 
for a single value of b (in this case b and b*, respectively). Therefore, the theorem 
follows. □ 

4. Uniqueness of the limit cycle for b < 817/1000 

In this section we will prove the uniqueness of the limit cycle of system ([2]) when 
b < 0.817. The idea of the proof is to find a suitable rational Dulac function for 
applying the following generalization of Bendixson-Dulac criterion. 

Proposition 4.1. Consider the C^- differential system 

j x = P{x,y), , . 

\y = Q{x,y), ^'^^ 
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and let W C be an open region with boundary formed by finitely many algebraic 
curves. Assume that: 

(I) There exists a rational function V{x,y) such that 

^ ' dx^ ~^ dy^ ^ i^dx ~^ dy^ ^^^^ 

does not change sign on U. Moreover M only vanishes on points, or curves 

that are not invariant by the flow of f[T^ . 
(II) All the connected components ofhl\ {V = 0}, except perhaps one, say U, 

are simple connected. The component U, if exists, is 1-connected. 
Then the system has at most one limit cycle in U and when it exists is hyperbolic 
and it is contained in U. Moreover its stability is given by the sign of —VM on U. 

The above statement is a simplified version of tlie one given in [5] adapted to 
our interests. Similar results can be seen in [H [TOl [HI [25] . 



Remark 4.2. Looking at the proof of Proposition we also know that: 
(i) The Dulac function used in the proof is 1/V. 

(a) In the region U, the curve {V{x,y) = 0} is without contact for the flow 
of f[T^ . In particular, by the Bendixson-Poincare Theorem, the ovals of the 
set {V{x, y) = 0} must surround some of the critical points of the vector 
field. 

To give an idea of liow we fiave found the function V that we will use in our 
proof we will first study the van der Pol system and then the uniqueness in our 
system when b < 0.615. Although we will not use these two results, we believe 
that to start studying them helps to a better understanding of our approach. 

4.1. The van der Pol system. Consider the Van der Pol system 



X = y, 

y = —X + (6^ — 



,2 .2^„. (17) 



Due to the expression of the above family of differential equations, in order to 
apply Proposition 14. 1^ it is natural to start considering functions of the form 

V{x,y) = f2y^ + fi{x)y + fo{x). 

For this type of functions, the corresponding M is a polynomial of degree 2 in y, 
with coefficients being functions of x. In particular the coefficient of y"^ is 

f[{x)+f2{b'-x'). 

Taking fi{x) = (x^ — 3b'^)f2x/3 we get that it vanishes. Next, fixing /2 = 6, and 
imposing to the coefficient of y to be zero we obtain that fo{x) = 6x^ + c, for any 
constant c. Finally, taking c = 6^(36^ ^4), we arrive to 

Vb{x, y) = 6y^ + 2{x^ - 3b'^)xy + 6x'^ + 6^(36^ - 4). (18) 
From ([T6|) of Proposition 14.11 the corresponding M, which only depends on x, is 

Mb{x,y) = Ax" + fe'(36' - 4)(x' - 6'). 

It is easy to see that for b e (0,2/^3) ^ (0,1.15), Mb{x,y) > 0. Notice that 
Vb{x,y) = is quadratic in y and so is not difficult to see that it has at most 
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one oval, see Figure [8] for b = 1. Then we can apply Proposition 14.11 to prove the 
uniqueness and hyperbolicity of the limit cycle for these values of b. 




Figure 8. The algebraic curve Vb{x, y) = with 6=1. 

We remark that taking a more suitable polynomial Dulac function, it is possible 
to prove the uniqueness of the limit cycle for all values of b, see p. 105]. We 
have only included this explanation as a first step towards the construction of a 
suitable rational Dulac function for our system ([2]). 

4.2. System with b < 651/1000. By making some modifications to the func- 
tion Vb given by ffTS]) . we get an appropriate function for system ([2]). Consider 

Vb{x, y) =[2x^ + 66^(1 - b'^)x\y^ + 6(1 - b'^)y'^ + 2{x'^ - 3b'^)xy 
+ 6(1 -62)x2 + 62(3^2 _4)_ 

Computing the double discriminant A'^iVb) of the function Vb, introduced in Ap- 
pendix II, we get that 

A\Vb) = b\3b' - 4){b^ - l)i5(Pi9(62))2, 

where Pig is a polynomial of degree 19. By using for instance the Sturm method, 
we prove that the smallest positive root of A'^{Vb) is greater than 0.85. Therefore 
by Proposition 15.71 we know that for b G (0,0.85] the algebraic curve Vb{x,y) = 
has no singular points and therefore the set {Vb{x, y) = 0} G is a finite disjoint 
union of ovals and smooth curves diffeomorphic to open intervals. 
By applying Proposition 14. II to system ([2]) with V = Vb, we get that 

Mb{x, y) = 6[(2 - 3fe2)a;V - 26^(2 - b^)xY + (2 - b^)x^y^] + 2(2 - 3b^)x^ 
-36^(14 - 1562)a:2y2 + 1254(2 - b^)xy^ _ 52^4 - 9b'^)x^ 
+3b\2 - 3b^)y^ + b\A - 3b^). 

(19) 

In Subsection 15.51 of Appendix II we prove that Mb does not vanish on for 
b G (0,0.651]. Then by Remark 14.21 all the ovals of {Vb{x,y) = 0} must surround 
the origin, which is the unique critical point of the system. Since the straight line 
X = has at most two points on the algebraic curve Vb{x, y) = 0, it can have at 
most one closed oval surrounding the origin. Then by Proposition 14.11 it follows 
the uniqueness, stability and hyperbolicity of the limit cycle of system ([2]) for these 
values of the parameter b. 
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4.3. System 1^ with b < 817/1000. The hyperbola xy + 1 = will play an 
important role in the study of this case. We first prove a preliminary result. 

Lemma 4.3. Consider system (|2]). 

(I) For b 1 the hyperbola xy + 1 = is without contact for its flow. In 
particular its periodic orbits never cut it. 

(II) For 6=1 the hyperbola xy + 1 = Q is invariant for its flow and the system 
has not periodic orbits. 

Proof. Define y) = xy + 1 and set X = {P,Q) := {y, —x + — x'^){y + y^)). 
Simple computations give that for x 7^ 0, 

Therefore (I) follows and we have also proved that when b = 1, the hyperbola is 
invariant by the flow. 

(II) When 6 = 1, 

F.,P + FyQ = KF, (20) 

where K = K{x, y) = y"^ — x"^ — xy{xy — 1) is the so called cof actor of the invariant 
curve F = 0. 

Let us prove that the system has no limit cycle. Recall that the origin is repeller. 
Therefore if we prove that any periodic orbit F of the system is also repeller we 
will have proved that there is no limit cycle. 

This will follow if we show that 

div(X)(7(t))dt > 0, (21) 



where 7(t) := {x(t),y(t)) is the time parametrization of F and T = T(F) its period. 

To prove (12T!) notice that the divergence of X can be written as div(X) = 
3K + 2x'^ + l- 3xy. Then, 

T rT pT 

div(X)(7(t))dt = 3 / K{x{t),y{t))dt+ / {2x{tf + l)dt - 3 / x{t)y{t)dt. 



Observe that from (120]) we have that 

Kixit)vmdt - r F^i^it)Mt))x + FMt)Mt))y ,^ 

K[x[t),y[t))at - Fixit),yit)) ^ 

\n\F{x{t),y{t))\dt = \Yi\F{x{t),y{t))\ ^ 



and that 



Therefore 



dt 



x{t)y{t)dt= / x{t)x{t)dt 



x\t) 



2 

T rT 



T 



0. 



div(X)(7(t))cit = / {2x{tf + l)dt>Q, 
Jo 

as we wanted to see. □ 

Theorem 4.4. System forb G (0, 0.817] has at most one limit cycle. Moreover 
when it exists it is hyperbolic and attractor. 
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Proof. Based on the function Vb used in the Subsection l4.2l we consider the function 

Vb{x,y) = Vb{x,y) /{5 + Gb^^x"^), where 

Vb{x, y) = \ b'^x' + \ fei W + (1 + I ^'y' + (1 + 1 ^'2/ 
- (I + f + 2 ^16^ ^2y2 + (3 52 _ 3 ^4 ^ |1 ^6) 

+ (3_352^264)a;2-62(3_i54)a;^ + (3_352 + 264)y2 ^^^^ 

+§6'^ -2 62. 

We have added the non-vanishing denominator to increase a httle bit the range 
of values for which Proposition 14.11 works. Indeed, it can be seen that the above 
function, but without the denominator, is good for showing that the system has 
at most one hmit cycle for b < 0.811. 

To study the algebraic curve Vb{x,y) = we proceed like in the previous sub- 
section. The double discriminant introduced in Appendix II is 

A\Vb) = b'^'^Sb^ - 4)(46=^^ + 276^4 + iQgb^^ + 108)(Pi52(6'))', 

where P152 is a polynomial of degree 152. It can be seen that the smallest positive 
root of A2(Vfe) is greater than 0.88. Therefore by Proposition 15 . 71 we know that for 
b G (0, 0.88] this algebraic curve has no singular points. Hence the set {Vb{x, y) = 
0} C is a finite disjoint union of ovals and smooth curves diffeomorphic to open 
intervals. 

The function that we have to study in order to apply Proposition 14.11 is 

N Nb{x,y) 
^^(^^^) = 30(6W + 5)2 (23) 

where Nb{x, y) is given in ( 133|) of Subsection 15.61 The denominator of Mb is positive 
for all {x,y) e M^. By Lemma 14.31 we know that the limit cycles of the system 
must lay in the open region VL = n{a;?/ + 1 > 0}. In Subsection 15. 6l of Appendix 
II we will prove that A*";, does not change sign on the region Q and if it vanishes it 
is only at some isolated points. 

Notice also that the set {Vb{x,y) = 0} cuts the y-axis at most in two points, 
therefore by the previous results and arguing as in Subsection 14.21 we know that 
it has at most one oval and that when it exists it must surround the origin. 

Therefore we are under the hypotheses of Proposition 14.11 taking U = Q, and 
the uniqueness and hyperbolicity of the limit cycle follows. □ 

5. Phase portraits for b < 79/100 and b = 817/1000 

This section is devoted to find the relative position of the separatrices of the 
infinite critical points when b < 0.79 and when b = 0.817. The main tool will 
be the construction of algebraic curves that are without contact for the fiow of 
system ([2]). These curves are essentially obtained by using the functions (piix) := 
(j)i{x — b) / {x — bY and ipi{x) := (pi{l/x) where 0j and ipi are the approximations of 
order i of the separatrices of the infinite critical points, given in the expressions ([3]) 
and (jl]) of Theorem 12. H respectively. That is, we use algebraic approximations of 
Si and (Sj', for i = 1,2, 3. 

As usual for knowing when a vector field X is without contact with a curve of 
the form y = ijj{x) we have to control the sign of 

N^x) := {\/iy-^ix)),X)\^^^^^y 
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In this section we will repeatedly compute this function when ip{x) is either (pi{x), 
4>i{x) or modifications of these functions. 
We prove the following result. 

Proposition 5.1. Consider system Then: 

(I) For b < 79/100 the configuration of its separatrices is like (i) in Figure\^ 
Moreover it has an odd number of limit cycles, taking into account their 
multiplicities. 

(II) Forb = 817/1000 the configuration of its separatrices is like (v) in Figurel^ 
Moreover it has an even number of limit cycles, taking into account their 
multiplicities. 



Proof. (I) Consider the two functions 

1 

X 



^i{x) 



and 



1 - 1: 

X x^ 



which are the corresponding expressions in the plane {x, y) of the first and second 
approximation of the separatrix Si. 

If 6 < 1 then [ipi — ip2){x) = {b"^ — l)/x^ > for a; < 0. This implies that the 
separatrix Si in the (x, ?/)-plane and close to — oo is below the graphic of (pi{x). 
Moreover 

M^{x) = _(^!±lK^!lli) < for a; < 0. 

X'^ 

This inequality implies that the separatrix Si in the plane (x, y) cannot intersect 
the graphic of (pi{x) for x < 0, see Figure O 



y 














X 



Figure 9. Behavior of 5i for 6 < 1. 



Now, we consider the third approximation to the separatrices ^2 and S-^, that 
is we consider the first three terms in ([3]) . It is given by the graph of the function 

^ , , (x^ - 5te + 1362) 

'^^^^^ = 9b^ix-br ■ 



Let us prove that when 6 G (0, 1^2/3), the graphs of (pi{x) and 03 (x) intersect 
at a unique point, (xo,?/o) with xq < and yo > 0. For this is sufficient to show 
that the function {ifi — 03)(x) has a unique zero at some xq < 0. 

It is clear that \im^_^Q-{(pi — 03)(x) = +oo and we have that {ipi — 03) (—26) = 
(36^ - 2)/6fe3, then for b < {<^i - 03) (-26) < hence {<^i - 03) (x) has a 
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zero at a point xq with —2b < Xq < 0. Moreover this zero is unique because the 
numerator of [ipi — (pzli^) is a monotonous function. 

It also holds that V(?/-</)3(x)) = (-03(x), 1) where = {-x + 7b) /^b'^{-x + 

6)'^ is a positive function for x < 0, and a simply computation shows that 

^^<^3(^) = 37296^3^ [(81^' + !)^' + (7296^ -4056^-11)^^' 

-9(1626^ - 1086^ - 7)62a;2 + (7296^ + 4056^ - n%)b^x 
-13(816^ - 20)6^] . 

To control the sign of we compute the discriminant of its numerator with 
respect to x. It gives d\s{N^.j^{x),x) = fc^^P22(&^), where P22 is a polynomial of 
degree 22 with integer coefficients. 

By using the Sturm method we obtain that P22(&^) has exactly four real zeros. 
By Bolzano theorem the positive ones belong to the intervals (0.7904, 0.7905) and 
(2.6,2.7). 




Figure 10. Behavior of Si and ^3 for b < 0.79 

If we fix bo < 79/100 then bo < ^/2/3 and moreover according to previous 
paragraph the graphics of fi{x) and (psi^x) intersect at a unique point {xo,yo) 
with xo < and yo > 0. Furthermore, -^(^0) > in (xq, 60) and Nfj,^ < in 
(xo, b) for all b G (0, bo]- Therefore the vector field associated to (|2]) on these curves 
is the one showed in Figure fTOl( a). 

From Figure [TOl(a) it is clear that the separatrix Si cannot intersect the set 
fl = {(x, (y9i(x))| — 00 < X < xo} U {(x, 03(x))|xo < X < bo}- Moreover, since 
the separatrix S2 forms an hyperbolic sector together with 1S3 we obtain that Si 
cannot be asymptotic to the line x = bo- Hence we must have the situation showed 
in Figure [Tni(b). We know that the origin is a source and from the symmetry of 
system ([2]) we conclude that for b < 0.79 the system has an odd number of limit 
cycles (taking into account multiplicities) and the phase portrait is the one showed 
in Figure [TT] 

(II) We start proving the result when b = bo '-= 89/100 because the method that 
we use is the same that for studying the case b = 817/1000, but the computations 
are easier. Recall that we want to prove that the configuration of separatrices is 
like (v) in Figure [31 That the number of limit cycles must be even (taking into 
account multiplicities) is then a simply consequence of the Poincare-Bendixson 
Theorem, because the origin is a source. 
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Figure 11. For < 6 < 0.79, system has at least one limit 
cycle and phase portrait (i) of Figure [T] . 

We consider the approximation of eight order to ^2 and 1S3 given by the graph 
of the function (t>s{x). 

By using again the Sturm method it is easy to see that N^^{x) < for x G 
(607 Xq), where Xq = 1.924 is a left approximation to the root of the function (psix)^ 
and N^g{x) > for a; G (0:1,60)5 where Xi = —2.022 is a right approximation 
to the root of the function N^g{x). That is, we have the situation shown in 
Figure [121 (a). Now, we consider the function (t>s{x) = (psi^) — is clear 

that (08 — (p8){x) > 0. We have iV^^(x) > for x G {bQ,X2) where X2 = 1.6467 is 

a left approximation to the root of the function 08(2;), moreover the line x = X2 
is transversal to the vector field for y > 0, thus the separatrix 1S3 intersects the 
X-axis at a point x of the interval (x2,Xo), see again Figure [T2l (a). 




X 



(a) (b) 

Figure 12. Behaviour of ^2 and ^3 for b G {0.817, 0.89}. 

At this point, the idea is to show that S2 intersects the x-axis at a point x, with 
— X2 < X < 0. For proving this, we utilize the Fade approximants method, see [3]. 
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Recall that given a function f{x), its Pade approximant Pdin,m]{f){x,Xo) of 
order {n,m) at a point Xq, or simply Pd[n,m\{f){x) when xq = 0, is a rational 
function of the form Fn{x)/Gm{x), where F„ and Gm are polynomials of degrees 
n and m, respectively, and such that 

0((x-xo)"+'"+^). 



fix)- 



G„ 



[X 



Consider the Pade approximant Pd[3^3](08). It satisfies that Pd[3_3]((/)8)(0) = (psi^) 
and by the Sturm method it can be seen that there exists X3 < such that 
Pd[3^3]((/)8)(a;3) = 0, Pd[3^3](08) is positive and increasing on the interval (2:3, 0) and 
a left approximation to X3 is —1.595. Moreover it is easy to see that iVpdjg 3^(^s){x) > 
Oforx G (x3,0). Therefore 1S2 cannot intersect neither the graph of ?/ = Pd[3 3](08)(3;) 
in (x3, 0) nor the graph of 08(3;) in [0, bo). Hence S2 intersects the x-axis in a point 
X contained in the interval (X3, 0). This implies that —X2 < 5; < as we wanted to 
see, because —X2 < X3. Hence the behavior of the separatrices is like Figure [T^(b). 
See also Figure fT3l 




Figure 13. Behavior of 52,53,5^ and S'^ for b G {0.817,0.89}. 

When feo = 817/1000 we follow the same ideas. For this case we consider 
the functions (pie{x) and 4>ig{x) = 4>m{x) — 1/(96'^). Recall that the graphic of 
4>ig{x) is the sixteenth order approximation to ^2 and ^3. It is not difficult to 
prove that A'^^^^ > on the interval {bo,X2), with X2 = 1.6421 and since the line 
X = X2 is transversal to X for y > 0, ^3 intersects the x-axis at a point x > X2- 
Also we have that A^^^g > on the interval (— 3/100, feo) and using the Pade 
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approximant Pd[5^i](0i6)(x, —3/100) we obtain that ^2 intersect to the x-axis in a 
point X G (xs, 0) with > —1.638. This imphes that —X2 < x < as in the case 
b = 0.89. Hence we have the same situation that in Figure [T51 □ 

Remark 5.2. As it is shown in the proof of Theorem the values 0.79 and 
0.817, obtained in the previous proposition, provide a lower and an upper bound 
for b*. We have tried to shrink the interval where b* lies using higher order ap- 
proximations of the separatrices, but we have not been able to diminish its size. 



Appendix I: The Descartes method 

Given a real polynomial P{x) = a^x"^ + ■ ■ ■ + aix + and a real interval 
/ = (a, /3) such that P{a)P{(3) ^ 0, there are two well-known methods for knowing 
the number of real roots of P in /: the Descartes rule and the Sturm method. 

Theoretically, when all the G Q and «, /3 G Q, the Sturm approach solves 
completely the problem. If all the roots of P are simple it is possible to associate 
to it a sequence of n + 1 polynomials, the so called Sturm sequence, and knowing 
the signs of this sequence evaluated at a and /3 we obtain the exact number of real 
roots in the interval. If P has multiple roots it suffices to start with P/ (gcd(P, P')), 
see [231 Sec. 5.6]. 

Nevertheless when the rational numbers have big numerators and denominators 
and n is also big, the computers have not enough capacity to perform the compu- 
tations to get the Sturm sequence. On the other hand the Descartes rule is not 
so powerful but a careful use, in the spirit of bisection method, can many times 
solve the problem. 

To recall the Descartes rule we need to introduce some notation. Given an 
ordered list of real numbers [6o, • • • , bn-i, bn] we will say that it has C changes 
of sign if the following holds: denote by [cq, Ci, . . . , Cm-i, Cm], m < n the new list 
obtained from the previous one after removing the zeros and without changing 
the order of the remaining terms. Consider the m non-zero numbers 5i := qq+i, 
i = 0, . . . , m — 1. Then C is the number of negative 5i. 

Theorem 5.3 (Descartes rule). Let C be the number of changes of sign of the list 
of ordered numbers 

[Oo, Ctl, a2, . . . , Ctn-l, 0"n\- 

Then the number of positive zeros of the polynomial P{x) = a„a;" -|- ■ ■ ■ -|- aix + oq, 
counted with their multiplicities, is C — 2k, for some G N U {0}. 

Corollary 5.4. With the notations of Theorem 15. 31 if C = then P{x) has not 
positive roots and if C = 1 it has exactly one simple positive root. 

In order to apply Descartes rule to arbitrary open intervals we introduce the 
following definition: 

Definition 5.5. Given a real polynomial P{x) and a real interval [a, /3) we con- 
struct a new polynomial 

N^{p){x) ■.= {x+ir^^p (^^). 
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We will call N^{P), the normalized version of P with respect to {a, (3). Notice 
that the number of real roots of P{x) in the interval is equal to the number 

of real roots of N^[P)[x) in (0, oo). 

The method suggested in [12] consists in writing (a, /3) = IJi=i(Q^j! ctj+i)) with 
a = ai < a2 < ■ ■ ■ < Oik < ctfc+i = /3 in such a way that on each it 
is possible to apply Corollary 15.41 to the normalized version of the polynomial. 
Although there is no systematic way of searching a suitable decomposition, we 
will see that a careful use of these type of ideas has been good enough to study 
the number and localization of the roots for a huge polynomial of degree 965, see 
Subsection 15.61 in Appendix II. 

Appendix II: A method for controlling the sign 

OF POLYNOMIALS IN TWO VARIABLES 

The main result of this appendix is a new method for controlling the sign of 
families of polynomials with two variables. As a starting point we prove a simple 
result for one-parameter families of polynomials in one variable. 

Let Gb{x) be a one-parametric family of polynomials. As usual, we write Ax{P) 
to denote the discriminant of a polynomial P{x) = a^x" -!-■■■ + aix + Qq, that is, 

A,.(P) = (-l)'^lRes(P(x),P'(x)), 
an 

where Res(P, P') is the resultant of P and P'. 
Lemma 5.6. Let 

Gbix) = ^„,(6)x" + + ■ ■ ■ + giib)x + goib), 

be a family of real polynomials depending also polynomially on a real parameter b 
and set f2 = M. Suppose that there exists an open interval / C M such that: 

(i) There is some bo G /, such that G'5(,(x) > on Q. 

(ii) For all be I, A^{Gb) 7^ 0. 
(Hi) For all b e I, gn{b) ^ 0. 

Then for all b e I, Gb{x) > on Q. 

Moreover if Q = Qb = (c(&), 00) for some smooth function c{b), the same result 
holds changing Q by this new Qb if wg add the additional hypothesis 

(iv) For all be I, Gb{c{b)) ^ 0. 

Proof. The key point of the proof is that the roots (real and complex) of Gb depend 
continuously of b, because gn{b) 7^ 0. Notice that hypotheses (iii) and (iv) prevent 
that moving b some root enters in Q either from infinity or from the boundary 
of Q, respectively. On the other hand if moving b some real roots appear from C, 
they do appear trough a double real root that is detected by the vanishing of 
Ax{Gb)- Since by item (ii), Arc{Gb) 7^ no real root appears in this way. Hence, 
for all b e I, the number of real roots of any Gb is the same. Since by item (i) for 
b = bo, Gbo > on Q, the same holds for all b e L □ 

To state the corresponding result for families of polynomials with two variables 
inspired in the above lemma, see Proposition I5.12[ we need to prove some results 
about the iterated discriminants (to replace hypothesis (ii) of the lemma) and to 
recall how to study the infinity of planar curves (to replace hypothesis (iii)). 
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5.1. The double discriminant. Let F{x,y) be a complex polynomial on C^. 
We write F as 

F{x, y) = any"" + a„-i2/""^ + a„_2|/""^ + . . . + ai|/ + ao, (24) 
where = ai{x) G C[x]. Then 

A,(F) = (-l)^^Res(F,9F/9y), 

and this resultant can be computed as the determinant of the Sylvester matrix of 
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We will write A^.j.(F) = l\x[/\y{F)). Analogously we can compute A^,^(F). 
This so called double discriminant plays a special role in the characterization of 
singular curves of {F{x, y) = 0} and it is also used in applications, see for instance 
[H [121 122] • In particular we prove the following result. 

Proposition 5.7. Let F{x,y) be a complex polynomial on C^. If{F{x,y) = 0} C 
has a singular point, that is, if there exists a point (xo,?/o) ^ such that 
F{xo,yo) = dF{xo,yo)/dx = dF{xo,y^)/dy = 0, then Al^{F) = Al^{F) = 0. 

Proof. We write F{x,y) in the form flMj) . Without lost of generality we assume 
that {xq, yo) = (0, 0). Then from the assumptions it follows that ao(0) = ao(0) = 
and ai(0) = 0, that is, ao(a;) = x'^ao{x) and ai{x) = xai{x), with both also 
polynomials. 

By using the Sylvester matrix S defined above, we have that 

det S = (-l)"ao det(^(2n - 1 | n - 1)) + ai det(5(2n - 1 | 2n - 1)), (25) 

where S{i \ j) means the matrix obtained from S by removing the i-th row and 
the j-th column. 

Notice that the elements of the last row of S{2n — 1 | 2?t, — 1) are only 0, ao 
and ai. Therefore, developing the determinant of this matrix from this row we get 
that det(S'(2?7. — 1 | 2r2 — 1)) = xQ{x), for some polynomial Q{x). 

Hence, by using ( 125|) . we get that det S = x'^P{x) with P(x) another polynomial. 
This implies that Ay{F) has a double zero at x = and hence Ay^^{F) = 0. 

Analogously we can prove that Ax{F) has a double zero aX y = and hence 
A^,,(F) = 0. □ 

Corollary 5.8. Consider a one-parameter family of polynomials Ff,{x,y), de- 
pending also polynomially on b. The values of b such that the algebraic curve 
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Fh{x,y) = has some singular point in have to be zeros of the polynomial 

A\F,) :=gcd {AlyiF,),Al,iF,)). 

By simplicity we will also call the polynomial A^(F;,), double discriminant of 
the family Fb{x, y). As far as we know the above necessary condition for detecting 
algebraic curves with singular points is new. 

Remark 5.9. (i) Notice that if in Corollaru \5.^ instead of imposing that for h G /, 
A2(F;,) ^ 0, it suffices to check only that either Al y{Fb) or .^(F^) ^ 0. 

(a) The converse of the Proposition \5. 7| is not true. For instance if we consider 
the polynomial F{x, y) = x^y^ + x + 1 then Al^{F) = Aly{F) = 0, however 
Fx{x,y) = 3x^y^ + 1 and Fy{x,y) = 3x^y'^ hence {F{x,y) = 0} does not have 
singular points. 

(Hi) Sometimes Ay^{F) ^ Aly{F). For instance this is the case when F = 
y"^ + x^ + hx^ + hx because 



Ai y{F) = -1105926^(6 - 4)(6 - and A;_^.(F) = 2566^(6 - 4). 
Notice that A^{F) = b^{b-A). 
5.2. Algebraic curves at infinity. Let 

F{x, y) = F\x, y) + F\x,y) + --- + y) 
be a polynomial on of degree n. We denote by 

y, z) = z^F\x, y) + z''-^F\x, y) + ■ ■ ■ + F^ix, y) 

its homogenization in MP^. 

For studying -F(x, y, z) in MP^ we can use its expressions in the three canonical 
charts of MP^, {[x : y : 1]}, {[x : 1 : 2;]}, and {[1 : y : 2;]}, which can be identified 
with the real planes {(x, 2;)}, and {{ij.z)} respectively. Of course the 

expression in the chart {[x : y : 1]}, that is, in the (x,|/)-plane is precisely F{x,y). 

We denote by Fi{x,z) and F2{y,z) the expressions of the function F in the 
planes {{x,z)} and {{y,z)}, respectively. Therefore Fi{x,z) = F{x,l,z) and 
F2{y,z) = F{l,y,z). 

Let [x* : y* : z*] E MP^ be a point of {F = 0}. If z* ^ 0, then [x* : y* : z*] 
corresponds to a point in M^, otherwise it is said that [x* : : 0] is a point of 
F at infinity. Notice that the points at infinity of F correspond to the points 
[x* : y* : 0] where [x*, y*) 7^ (0, 0) is a solution of the homogeneous part of degree 
n of F, 

Hn{F{x,y))=F^{x,y), 

that is F"-{x*,y*) = 0. Equivalently, these are the zeros of Fi{x,0) and ^2(^,0). 
In other words, [x* : : 0] is a point at infinity of F if and only if x* /y* is a zero 
of Fi{x, 0) = 1) or y*/x* is a zero of ^2(7/, 0) = F"(l, y). 

Let r2 C be an unbounded open subset with boundary dfl formed by finitely 
many algebraic curves. It is clear that this subset can be extended to MP^ We 
will call the adherence of this extension fl. When a point at infinity of F is also 
in 0, for short we will say that is a point at infinite which is also in Q. 
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5.3. Isolated points of families of algebraic curves. To state our main re- 
sult we need explicit conditions to check when a point of a real algebraic curve 
G{x, y) = is isolated. Recall that it is said that a point p G on the curve is 
isolated if there exists an open neighborhood U of p, such that 

un{{:x,y)ew? G(x,y) = o} = p. 

Clearly isolated points are singular points of the curve. Next result provides an 
useful criterion to deal with this question. 

Lemma 5.10. Let G{x,y) be a real polynomial. Assume that (0,0) G {G{x,y) = 
0} and that there are natural numbers p,q and m, with gcd(p, g) = 1, and a 
polynomial G° satisfying G'^{e'''X, e'^Y) = Y), and such that for all e > 0, 

Gie^X, e'^Y) = Y) + e™+^G^(X, Y, e), 

for some polynomial function G^ . If the only real solution of G^{X,Y) = is 
{X,Y) = (0,0), then the origin is an isolated point of G{x,y) = 0. 

Proof. Assume without loss of generality that G° > 0. We start proving that 
K := {{x,y) G : G^{x,y) = 1} is a compact set. Clearly it is closed, so it 
suffices to prove that it is bounded. Since G^ is a quasi- homogeneous polynomial 
we know that there exists a natural number mo such that m = rriQpq and G^{x, y) = 
Pmo{x'^,y^), where P^o is a real homogeneous polynomial of degree mo. The fact 
that the only real solution of the equation G^{x,y) = is x = y = implies that 
Pmo has not linear factors when we decompose it as a product of real irreducible 
factors. Hence mo is even and Pmg{x,y) = YYh^i'^ i^i^'^ + BiXy + G^y"^), with 
Bf — iAiGi < 0. As a consequence, 

mo/ 2 

4 = 1 

Assume, to arrive to a contradiction, that K is unbounded. Therefore it should 
exist a sequence {{xn, yn)}, tending to infinity, and such that G^{xn, y-n) = 1- But 
this is impossible because the conditions Bf — AAiGi < 0, z = 1, . . . ,mo/2, imply 
that all the terms Aix'^'' + BiX'^y^ + G^y"^ in fl2B]) go to infinity. So K is compact. 

Let us prove that (0,0) is an isolated point of {{x,y) G : G{x,y) = 0}. 
Assume, to arrive to a contradiction, that it is not. Therefore there exists a 
sequence of points {(x„, ?/„)}, tending to and such that G{xn,yn) = for all 
n G N. Consider G^{xn,yn) ='■ {gn)'^ > 0. It is clear that lim„^oo(5'n)'" = 0. Write 

{Xn,yn) = {{gnYUn, {gnYVn)- NoticC that 

[gnT = G'°(x„,2/„,) = G°((7>„,,^7>„) = [gnTG^Un^Vn). 

Then G'''(u„,i;„) = 1 and (u„,f„) G if, for all n G N. Therefore, taking a subse- 
quence if necessary, we can assume that 

Yim{un.Vn) = {u\v*) eK. (27) 

n— ^-oo 

We have that = G(x„,2/„,) = {gnT + {9nr^'G\un,Vn, gn). Dividing by {gn)"" 
we obtain that 1 = -|- gnG^{un, Vn, gn), and passing to the limit we get that 1 = 
which gives the desired contradiction. 

Notice that to prove that lim^^oo 9nG^{un, Vn, gn) = we need to know that the 
sequence {{un,Vn)} remains bounded and this fact is a consequence of flTTI) . □ 
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We remark that the suitable values p, q and m and the function G° appearing 
in the statement of Lemma 15.101 are usually found by using the Newton diagram 
associated to G. 

We also need to introduce a new related concept for families of curves. Consider 
a one-parameter family of algebraic curves Gb{x,y) = 0, 6 G /, also depending 
polynomially of b. Let (xo,|/o) ^ be an isolated point of Gb{x,y) = for all 
6 e /, we will say that {xo,yo) is uniformly isolated for the family Gb{x,y) = 0, 
6 G / if for each b ^ I there exist neighborhoods V C / and W C M^, of 6 and 
{xo,yo) respectively, such that for all 6 G V, 

{(x, G : Gb{x, y) = 0} n W = (xo, yo). (28) 

Next example shows a one-parameter family of curves that has the origin isolated 
for all 6 G M but it is not uniformly isolated for 6 G /, with G /, 

Gb{x, y) = {x' + y'){x' + y^ - b'){x - 1). (29) 

It is clear that the origin is an isolated point of {Gb{x, y) = 0} for all 6 G M, but 
there is no open neighborhood W of (0,0), such that fl28|) holds for any 6 in a 
neighborhood of 6 = 0. 

Next result is a version of Lemma 15.101 for one-parameter families. In its proof 
we will use some periodic functions introduced by Lyapunov in his study of the 
stability of degenerate critical points, see [U]. Let us recall them. 

Let u{ip) = Cs{(p) and v{ip) = Sn{{p) be the solutions of the Cauchy problem: 



u' = -v-'^ \ V' = u^''~\ n(0) = VVl' and v{0) = 0, 



v' 

where the prime denotes the derivative with respect to (/?. 

Then x = Cs{(p) and y = Sn((/?) parameterize the algebraic curve px'^'^+qy'^^ = 1, 
that is pCs^'^(9?) -|- gSn^^((y9) = 1, and both functions are smooth Tp^^-periodic 
functions, where 



T = Tp. = 2p-i/2</g-i/2p. 



2pJ ^ \ 2g 



^ + ^ 

2p ^ 2q 



and r denotes the Gamma function. 



Proposition 5.11. Let Gb{x,y) be a family of real polynomials which also de- 
pends polynomially on b. Assume that (0,0) G {Gb{x,y) = 0} and that there are 
natural numbers p,q and m, with gcd(p, g) = 1, and a polynomial satisfying 
GliePX, e'^Y) = e^^G^iX, Y), and such that for all e > 0, 

Gb{ePX,e''Y) = e'^G'',{X,Y) + e'^+'Gl{X,Y,e), 

for some polynomial function Gl . If for a// 6 G / C M, the only real solution of 
G°(X, y) = is {X,Y) = (0,0), then the origin is an uniformly isolated point of 
Gb{x, y) = for all b e I. 

Proof. Assume without loss of generality that G^ > 0. Let us write the function 
Gb{x,y) using the so-called generalized polar coordinates, 

x = //Cs((^), ?/ = p«Sn((^), for p G M+. 
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Then 

G,(x,y)=G,(/Cs(^),p'^Sn(^)) 

= p'"G°(Cs(y.),Sn(y,)) +p-+iG,i(Cs(y.),Sn(y,),p). (30) 

Using the same notation that in the proof of Lemma 15.101 with the obvious mod- 
ifications, we know from fl26|) that 

mo/2 

G'°(Cs((/.), Sn(y.)) = n {A{b) Cs^'^iv) + B,ib) Cs%^) Sn^ (y.) + ^(5) Sn'^iip)), 

i=l 

with all Bf{b) - 4Ai{b)Ci{b) < 0. Therefore, it is not difficult to prove that there 
exists two positive continuous functions, L{b) and U{b) such that 

0<L(6)<G°(Cs(^),Sn(^))<t/(6), 

due to the periodicity of the Lyapunov functions and the discriminant conditions. 
Dividing the expression (15U]) by we obtain that the points of {Gb{x,y) = 
(0,0)} \ {(0,0} are given by 

Gt{Cs{^), Sn((p)) +pGl{Cs{^), Sn(y,), p) = 0. (31) 

Fix a compact neighborhood of b, say V C /. Set L = min^-gy -^(^)- Then there 
exists 6 > such that for any ||(x,i/)|| < 6 and any 6 e V, 

|pG,i(Cs(^),Sn(^),p)|<L/2. 

Therefore (13T|) never holds in this region and 

{{x,y)eR' : G,{x,y) = 0}n{{x,y)eR^ : \\{x,y)\\ < 6} = (0,0), 
for all 6 G V, as we wanted to prove. □ 

Notice that, the fact that for all 6 G M, the origin of (129|) is isolated simply 
follows plotting the zero level set of Gb. Alternatively, we can apply Lemma [5.101 
with p = l,q = 1 and m = 2 to prove that the origin is isolated when b ^ and 
with p = q = 1 and m = 4 when 6 = 0. In any case. Proposition 15.111 can not be 
used. 

5.4. The method for controlling the sign. 

Proposition 5.12. Let Ff,{x,y) be a family of real polynomials depending also 
polynomially on a real parameter b and let Q G M."^ be an open connected subset 
having a boundary dQ formed by finitely many algebraic curves. Suppose that there 
exists an open interval / C M such that: 

(i) For some bo G /, Fb^{x,y) > on Q G M^. 
(ii) For all be I, A\Fb) ^ 0. 

(Hi) For all b E I, all points of Fb = at infinity which are also in Q do not 

depend on b and are uniformly isolated, 
(iv) For all be I, {Fb = 0}ndn = $. 

Then for all b G /, Fb{x, y) > on Q. 
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Proof. Consider the following set 

J ■.= {bel : Fb{x, y)>0 for all (x, y) G Q}. 

By hypothesis (i), J because G J. Consider now h = sup J. We want to 
prove that b G dL If this is true, arguing similarly with inf J the result will follow. 
We will prove the result by contradiction. So assume that b & L 
Notice that if Fi{x,y) takes positive and negative values on Q, by continuity 
this would happen for any b near enough to b. This is in contradiction with the 
fact that b is the supremum of J. Therefore, either Fi{x, y) > or Fi{x, y) > 
in Q. 

In the first case it is clear that a point {xq, yo) where Fi{xo, yo) = has to be a 
singular point of the curve {Fi{x, y) = 0}. Therefore, by Corollary 15. 8 [ A''{Fi) = 
which is in contradiction with (ii). 

In the second case it should exist a sequence of real numbers {6„}, with 6„ I b, 
and a sequence of points {(a:„, G fl such that lim„_).oo F^^^Xn, yn) = 0. 

If the sequence is bounded, renaming it if necessary, we arrive to a convergent 
sequence. Call (x, |/) G Vt its limit, where Vt denotes the adherence of Vt. Then 
Fiix.y) = 0. By hypothesis (iv), the point {x,y) ^ dfl and we also know that 
Fi{x,y) > on r2. Therefore we have a contradiction and the sequence {{xn,yn)} 
must be unbounded. 

This unbounded sequence can be considered in the projective space MP^ Then 
this sequence must converge to a point p of Fi{x,y) = at infinity, which is 
also in U. Since by hypothesis (iii) this point is uniformly isolated, there exists 
a neighborhood V of 6 and an open neighborhood W of p such that this point 
is the only real point in MP^ of the homogenization of Fb{x, y) = 0. This is in 
contradiction with the fact Fb„ {xn, yn) = for all n, and the result follows. □ 

5.5. Control of the sign of (|T9|) . In this subsection we will prove by using 
Proposition I5.12( that for b G (0, 0.6512), the function given in ( IT9|) is positive 
onn = M2. 

To check hypothesis (i), we prove that Ml/2 > for all M^. por this value, 

A/f , — 15 ^4„,2 _ 21 3„,3 , 21 2„,4 _ 123 2„,2 , 21 3 i 5 4 _ X ^.2 i 15 „,2 , 13 

We think M1/2 as a polynomial in x and y as a parameter and we apply Lemma 
li y = then M1/2 reduces to the polynomial {5/2)x' — (7/16)x^ + 13/64 which 
is positive on M. Now, we compute Ax{Mi/2) and we obtain a polynomial in the 
variable y of degree 20. By using the Sturm method it is easy to see that it does 
not have real roots. Moreover, the coefficient of is 5(3?/^ + l)/2 > 0. Therefore, 
Ml/2 > on M^, as we wanted to see. 

To check hypothesis (ii) we compute the double discriminant of and we 
obtain that A^. ^(Mf,) is a polynomial in b of degree 1028, of the following form 

Aly{Mb) =5320^^2 _ 2)40(3^2 _ 2)5^3^2 _ 4^^256 _ 454 _ 3^2 ^2)x 

X (66 _ 2b' - 3b' + 2){P2{b'))'{Pe{b'))\Ps2{b'))\Ps^{b'))^ 

where Pi are polynomials of degree i with rational coefficients. By using the Sturm 
method we localize the real roots of each factor of A'y{Mb) and we obtain that 
in the interval (0, 0.6512) none of them has real roots. In fact P32{b') has a root 
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in (0.6513, 0.6514) and that is the reason for which we can not increase more the 
value of b. Therefore Al y{Mb) ^ for all b e (0, 0.6512). 

Finally we have to check hypothesis (iii). Notice that in this case dQ = and 
so (iv) follows directly. 

The zeros at infinity are given by the directions 

neiMb) = QxYli^ - 3fc')x2 - 2b\2 - b^)xy + (2 - b^)y^] = 0. 

For \b\ < 0.7275 it has only the non-trivial solutions x = and y = 0. The 
homogenization of Mf, is 

Mb = 6[(2-3fe2)xV-262(2-&')xV + (2-fe')xV]+2(2-362)xV 
-3fe2(14 - 15fe2)a;2|/2^2 + i2b\2 - b^)xy'^z^ - 6^(4 - 9b^)x''z^ 
+3b\2 - 3b^)y^z^ + b\4: - 3b^)z'^, 

(32) 

and hypothesis (iii) is equivalent to prove that (0, 0) is an uniformly isolated 
singularity for Ml{x, z) = Mb{x, 1, z) and that (0, 0) is also an uniformly isolated 
singularity for M^^y, z) = Mf,{l, y, z). 

First we prove this result for Ml{x,z). From (132|) . 

Ml{x, z) = 6[(2 - 362)a;4 - 2b\2 - b^)x^ + (2 - b^)x^] + 2(2 - 3b^)x^z^ 
-3b\U - 15b^)x^z'^ + 126^(2 - b'^)xz'^ - b\4 - 9b'^)x^z^ 
+3b\2 - 3b^)z^ + 6^(4 - 3b'')z^. 

Hence, 

M^{e^X,eZ) = ^6(2 - b^)X^ + 126-^(2 - b^)XZ^ + 36^(2 - 3b^)Z^y^ + Oie'). 

The discriminant with respect to X of the homogeneous polynomial T{X, W) := 
6(2 - 62)X2 + 126^(2 - b'^)XW + 36^(2 - 3b'^)W^, where W = Z^, is 

Ax(T) = 72W%\b^ - 2){2b^ - Ab^ - 3b^ + 2). 

Since its smallest positive root is greater than 0.673 it holds for b G (0, 673) that 
T{X,W) = if and only if {X,W) = (0,0). Therefore by Proposition EIII] the 
point (0,0) is an uniformly isolated point of the curve M^{x,z) = 0, for these 
values of b. 

For the other point, since 

M^{y^ z) = 6[(2 - b^)y^ - 2b\2 - b^)y^ + (2 - 362)^2] + 2(2 - 362)^2 

-3fe2(14 - 1562)^2^2 ^ ^26^(2 - 62)^3^2 _ ^2^4 _ 9^2)^4 
+3b\2 - 362)^2^4 + 54^4 _ 3^2)^6^ 

we have that 

M^ieY, eZ) = 2(2 - 36^) [3Y^ + Z'^'^e^ + 0{e^), 

and the result follows for b G (0, ^2/3) ~ (0,0.816), by applying again the same 
proposition. 

So, we have shown that for b G (0, 0.6512) all the hypotheses of the Proposi- 
tion [5]T2] hold. Therefore we have proved that for b G (0,0.651], Mb{x,y) > for 
all {x,y) G M^. 



BIFURCATION VALUES 29 

5.6. Control of the sign of fl23l) . The numerator of the function Mi, given in 
f l2^ is a polynomial of the following form 

Nb{x, y) = /o(x, h) + 6)y + ^(x, + /3(x, 6)y3 + f^^^^ b)y\ (33) 

where 

fo{x, b) = 906^^^° - 15618(6620 - 5)x^ + 1561^(246^ - 596^ + 2A)x'^ 
-(37862^^ - 810622 + 360620 - 3006^ + 67 562 _ soo)^^ 
-1562(18622 - 24620 + 216^ - 4562 + 20)^2 - 7564(-4 + 362), 

/i(x, 6) = 180636x'^ + 12618(60616 + 506^^ + 186i° + 25)x^ - 206iO(366i2 
-54610 + 546^ - 3066 _ 2554 _ 9)^3 _ i80620(362 - 4)x, 

/2(x,6) = 270636x10 - 45618(6620 + 2618 -5)x8 + 3618(30620 + 120616 
+10061^1 - 90612 _^ 35^10 _^ 35054 _ 5^^552 _^ 335)^:6 _ (3606^6 
+3006^4 + 1086^0 + 22 14624 _ 3590622 _^ 3435520 _^ 350518 
-300616 - 25061"^ + 225612 - 906io - 9006^ + 135062 - 900)x4 
-62(468622 - 540620 - IO8O6I8 + 3006i6 + 25061^ + 906i° 
+18456^ - 307562 _^ 2475)x2 - 906"^ (462 - 5), 

/3(x, 6) = -180620(510 - 3)x7 + 3062 (663^ ^ g^so _ 24522 ^ ^8^20 _ 725I8 
-5610 + 15)x5 + 3062(24624 - 36622 ^ 72620 + IO6I6 + 56i2 
-2064 + 1562 - 60)x3 - 2064(36618 - 546i6 + 546i4 + 3O612 
+2561° + 966 _ 3054 _^ 4552 _ 90)3,^ 

U{x, 6) = 906^6x8 - 3618(30620 + 12O6I6 + 1006i4 + 366io - 25)x6 
+610(36 0626 + 30 0624 ^ ^93520 ^ 350512 _ 5^5^10 ^ 72O68 
-30066 _ 25054 _ 90)x4 + (-738624 ^ 108O622 - 1080620 
+300618 + 250616 + 315612 + 30064 _ 45052 _^ 900)^2 + 1566. 

We will prove that > on n := {{x,y) : xy + 1 > 0} for all b e (0, 0.817] 
and if it vanishes this only happens at some isolated points. We will use again 
Proposition 15.121 Notice that dfl = {{x,y) : xy + 1 = 0}. 




Figure 14. Curves A^^ = and xy + 1 = with b = 0.817. 

It is not difficult to verify that {Nb{x, y) = 0} H {xy + 1 = 0} = for 6 G 
(0, 0.8171), see Figure IT^ It suffices to see that for these values of b, and x ^ 0, the 
one variable function Ni,{x,l/x), never vanishes. We skip the details. Therefore 
hypothesis (iv) is satisfied. 
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For proving that hypothesis (ii) of Proposition [5TT2] holds we compute the double 
discriminant Ay^{Nf,). It is an even polynomial in b, of degree 21852, of the 
following form 

(34) 

where Pi are polynomials of degree i with rational coefficients. By using the Sturm 
method it is easy to see that its first 4 factors do not have real roots in (0, 0.8171). 
We replace 6^ = t in the next three polynomials to reduce their degrees and 
we obtain Vi{t) := Psseit), V2{t) := P^siit), and V-sit) := Pg65{t). It suffices to 
study their number of real roots in (0,0.6678], because 0.6678 > (0.8171)^. Our 
computers have no enough capacity to get their Sturm sequences. Therefore we 
will use the Descartes approach as it is explained in Appendix I. 

We consider first the polynomial Vi{t). Its normalized version Nq'^^{Vi) has 
all their coefficients positive. Therefore Vi{t) has no real roots in (0,0.68) as we 
wanted to see. 

Applying the Descartes rule to the normalized versions of V2{t), Nq'^^^{V2), 
N^-Ml{V2) and A^o.lel (^2), we obtain that the number of zeros in the intervals 
(0, 0.561), (0.561, 0.811) and (0.562, 0.812) is 0, 1 and respectively. That is, there 
is only one root of V2(t) in (0, 0.812), it is simple and it belongs to (0.561, 0.562). 
Refining this interval with Bolzano Theorem we prove that the root is in the 
interval (0.5617,0.5618). 

Finally to study V^it) we consider N^^^^^Vs), iVn/2o(p3) and NyH^Vs). By 
Descartes rule we obtain that the number of zeros of P3 in the corresponding 
intervals is 0, 1 and 1 or 3, respectively. By Bolzano Theorem we can localize 
more precisely these zeros and prove that in the last interval there are exactly 3 
zeros. So we have proved that the polynomial V3 has exactly 4 zeros in the interval 
(0, 52/75) ~ (0, 0.693), and each one of them is contained in one of the following 
intervals 

(0.5614,0.5615), (0.6678,0.6679), (0.6690,0.6700), (0.6870,0.6880). 

In brief, for t G (0,0.6678] the double discriminant Ay^x{Nb) only vanishes at two 
points t = h and t = with h e (0.5614, 0.5615) and t2 e (0.5617,0.5618). 
Therefore we are under the hypothesis (ii) of Proposition 15.121 for b belonging to 
each of the intervals (0,6i), (fei,&2) and (62, 0.8171), where 

b^.-^r^ 0.749301, 62 := ~ 0.749478. 

To ensure that on each interval we are under the hypotheses (i) of the proposition 
we prove that does not vanish on Q for one value of b in each of the above 
three intervals. We take 

1 7494 3 

2^(0,61), (61,62), and -€(62,0.8171). 
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We study with detail the case b = 1/2. The other two cases can be treated similarly 
and we skip the details. So we have to study on Q the sign of the function 

/V , ^ 135 ^10. .2 I 45 8„,4 , 45 10 , 117964485 8„,2 

1/2 34359738368 " ^ 34359738368 " ^ 34359738368 ^ 137438953472 f 

I 138195 7 3 I 39253779 6 4 , 39321555 8 , 45 7 

~'~ 268435456 tl "r 137438953472-^ " 137438953472-^ 17179869184"^ » 

I 320504301 6„,2 , 1932072223485 5,, 3 _ 906074381 4 4 , 645 6 
137438953472-'' W 17179869184 » 8589934592"'' » 1048576-^ 

I 1229859 5„, , 5315442024413 4 2 _ 1808748465 3,, 3 , 6763995071 2 4 
1073741824 -'' » 8589934592 -'' f 4194304 U ~^ 8388608 » 

I 1258289751 4 , 55625 3 _ 1910154937 2,, 2 , 26361865 3 , 15 4 
8388608 262144"'' f 4194304 U ~<~ 262144 "''» "■"64*' 

316538295 2 , 585 , 45 2 , 975 

8388608 1048576 -''w ~'~ 2 » 64 ' 

We consider N1/2 as a polynomial in x with coefficients in 'R[y] and we apply 
Lemma I^TBl with fij^ = (— 00) when y > and f2o = (—00,00). Notice 
that for the symmetry of the function there is no need to study the zone y < 
because Ni/2{—x, —y) = Ni/2{x,y). We introduce the following notation Sy{x) : = 
Ni/2{x,y). We prove the following facts: 

(i) If we write Sy{x) = J2i=i^^i{y)^\ then sio(y) = A;(l + 3y^) for some k G 
Therefore sio{y) > for all e M. 

(ii) li y = then So{x) is an even polynomial of degree 10 and it is easy to see 
that So{x) > over M. 

(iii) We already know that {Sy{x) = 0} n 9fi = 0. 

(iv) Some computations give that 

A^iSy) = Ps,{y'), 

where P35 is a polynomial of degree 35. Moreover, using once more the Sturm 
method, we get that P-i^iy"^) has only two positive roots Q < yi < y2, with 
7/1 ~ 0.588423 and y2 ~ 6065.2946. From this result it is easy to prove that: 

(a) If ?/ e [0,yi) U (?/2, 00), then Sy{x) > 0. 

(b) If ?/ e (2/1,2/2), then Sy{x) has only two real roots, say xi{y) < X2{y), 
and none of them belongs to the interval {—1/y, 00). So Sy{x) > on 
(-l/y,oo). 

(c) If y G {7/1,2/2}, then Sy{x) has only a real root, Xi{y), which is a double 
root and Xi{y) ^ ( — l/y, 00). So, again Sy{x) > on (— l/y, 00). 

Thus, by Lemma the function N1/2 is positive on {x,y) G fl, as we wanted 
to see. In fact, its level curves are like the ones showed in Figure [HJ The straight 
lines y = yi and 2/ = 2/2 correspond to the lower and upper tangents to the oval 
contained in the second quadrant. 

To be under all the hypotheses of Proposition 15.121 it only remains to study the 
function Nf, at infinity. We denote by iVb(a;, y, z) its homogenization in MP^ and by 
Nl{x, z) and N^{y, z) the expressions of the function A'';, in the planes {(x, z)} and 
{(2/, 2;)}, respectively. Since T-LuiNb) = 90b^^x^y'^[3x^ + 2/^], the only non-trivial 
solutions of 'Hi2{Nb) = are x = and y = 0. Hence these directions give rise to 
two points of Ni) at infinity which are also on the region Q. They correspond to 
the points (0, 0) of the algebraic curves Nl{x, z) = and N^{y, z) = 0. We have 
to prove that both points are uniformly isolated. 
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Similarly that in the previous subsection, we write 
Nl{eX,eZ) = {%Qb^^X^ - 36^^(306^° + 1206^^ + 1006" + 366^" - 2b)X^Z'^ 

+ 6i"(360626 + 3006^4 + 19862° + 3606^^ _ g^g^io ^ ^206^ - 3006^ - 2506"* - 9{))X'^Z^ 
+ (-73852^ + 10806^2 - 10806^° + 3006^*^ + 2506^*^ + 3156^^ ^ ^QQ^i _ 45952 ^ qqq)x'^Z^ 
+ 1566z**)e« + 0(e9) 

and _ 

Nl{eY,eZ) = m^\3Y^ + Z^)e^ + 0{e'). 

By Proposition I5.1H for the second algebraic curve it is clear that for all 6 > 
(0, 0) is an isolated point. 

For studying the first one we denote by R{X, Z) the homogenous polynomial 
accompanying and we obtain that 

Ax{R{X,Z)) = Z'%''%Pnm\ 

for some polynomial P71 of degree 71 and integer coefficients. Since the smallest 
positive root of this polynomial is greater that 0.92 we can easily prove that for 
b < 0.92, R{X, Z) = if and only if X = Z = 0. Therefore we can use again 
Proposition 15.111 and prove that (0, 0) is an uniformly isolated point of the curve 
for these values of b. 
So, if we write 

(0, 0.8171) = (0, 61) U {61} U (61, 62) U {62} U (62, 0.8171), 

we can apply Proposition 15.121 to each one of the open intervals to prove that for 
b e (0, 0.817] \ {bi, 62} it holds that Nb{x, y) > for all {x, y) in VL. By continuity, 
for the two values b G {bi,b2}, we obtain that Nh{x,y) > 0. Since Ay^Nh) ^ 
either it is always positive or it vanishes only at some isolated points, as we wanted 
to prove. 

It can be seen that for b ~> b ^ 0.81722, Nh{x, y) changes sign on Q because there 
appears one oval in the set {Nh{x, y) = 0}. The value 6^ ~ 0.6678492 corresponds 
to the root of V3 in the interval (0.6678, 0.6679) that has appeared in the proof as 
a root of the double discriminant. 
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